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Abstract. We present a systematic theory of the phonon Hall effect in a ballistic 
crystal lattice system, and apply it on the kagome lattice which is ubiquitous in 
various real materials. By proposing a proper second quantization for the non-Hermite 
Hamiltonian in the polarization-vector space, we obtain a new heat current density 
operator with two separate contributions: the normal velocity responsible for the 
longitudinal phonon transport, and the anomalous velocity manifesting itself as the 
Hall effect of transverse phonon transport. As exemplified in kagome lattices, our 
theory predicts that the direction of Hall conductivity at low magnetic field can be 
reversed by tuning temperatures, which we hope can be verified by experiments in 
the future. Three phonon-Hall-conductivity singularities induced by phonon-band- 
topology change are discovered as well, which correspond to the degeneracies at three 
different symmetric center points, F, K, X, in the wave-vector space of the kagome 
lattice. 
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1. Introduction 

In recent years, phononics, the discipline of science and technology in processing 
information by phonons and controlling heat flow, becomes more and more exciting 
[H [2] . Various functional thermal devices such as thermal diode [3] , thermal transistor 
[1], thermal logic gates [5] and thermal memory [6], etc., have been proposed to 
manipulate and control phonons, the carrier of heat energy and information. And 
very recently, similar to the Hall effect of electrons, Strohm et al. observed the 
phonon Hall effect (PHE) - the appearance of a temperature difference in the direction 
perpendicular to both the applied magnetic field and the heat current flowing through 
an ionic paramagnetic dielectric sample [7] , which was confirmed later in Ref . [8j . Such 
observation of the PHE is really surprising because phonons as charge- free quasiparticles, 
different from electrons, cannot directly couple to the magnetic field through the Lorentz 
force. Since then, several theoretical explanations have been proposed [9l [TOl [TTl |12] to 
understand this novel phenomenon. From the work of the PHE in four-terminal nano- 
junctions and the phonon Hall conductivity in the two-dimensional periodic crystal 
lattice, we know that the PHE can exist even in the ballistic system. 

Geometric phase effects [131 El ^"^^ fundamentally important in understanding 
electrical transport property in quantum Hall effect [151 US], anomalous Hall effect [T71 
HE], and anomalous thermoelectric transport [19]. It is successful in characterizing 
the underlying mechanism of quantum spin Hall effect [201 EI]- Such an elegant 
connection between mathematics and physics provides a broad and deep understanding 
of basic material properties. Although there is a quite difference between phonons and 
electrons, we still can use the topological description to study the underlying properties 
of the phonon transport, such as topological phonon modes in dynamic instability of 
microtubules [22] and in filamentary structures [23], Berry-phase- induced heat pumping 
|24] . and the Berry-phase contribution of molecular vibrational instability [25] . 

The topological nature of the PHE is recently studied in Ref. [26] , where a general 
expression for phonon Hall conductivity is obtained in terms of the Berry curvatures 
of band structures. In Ref. [26], the authors find a phase transition in the PHE of the 
honeycomb lattice, explained from topological nature and dispersion relations. From the 
Green-Kubo formula and considering the contributions from all the phonon bands, the 
authors obtain the general formula for the phonon Hall conductivity. Then by looking at 
the phases of the polarization vectors of both the displacements and conjugate momenta 
as a function of the wave vector, a Berry curvature can be defined uniquely for each band. 
Combining the above two steps, at last the phonon Hall conductivity can be written in 
terms of Berry curvatures. Such derivation gives us a clear picture of the contribution 
to the phonon Hall current from all phonon branches, as well as the relation between the 
phonon Hall conductivity and the geometrical phase of the polarization vectors, which 
thus helps us to understand the topological picture of the PHE. However, the process 
of going from the Berry phase to the heat flux and the phonon Hall conductivity looks 
not very clear and natural. 
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We know that for the Hall effect of the electrons, in addition to the normal velocity 
from usual band dispersion contribution, the Berry curvature induces an anomalous 
velocity always transverse to the electric field, which gives rise to a Hall current, thus 
the Hall effect occurs [H] . For the magnon Hall effect [27] recently observed, the authors 
also found the anomalous velocity due to the Berry connection which is responsible for 
the thermal Hall conductivity. Therefore in this article we will derive the theory of 
the PHE in a more natural way where the Berry phase effect inducing the anomalous 
velocity contributes to the extra term of the heat current. Thus the Berry phase effect 
is straightforward to take the responsibility of the PHE. 

A kagome lattice, composed of interlaced triangles whose lattice points each have 
four neighboring points [29], becomes popular in the magnetic community because the 
unusual magnetic properties of many real magnetic materials are associated with those 
characteristic of the kagome lattice [30] . The schematic figure of kagome lattice is shown 
in Fig. [U In this paper we also apply the PHE theory to the kagome lattice to investigate 
whether the mechanism of the phase transition found in Ref. [26] is general and how 
the phonon Hall conductivity, Chern numbers and the dispersion relation behave and 
relate to each other. 

In this paper we organize as follows. In Sec. 121 we give a new systematic derivation of 
the theory of the PHE in terms of Berry curvatures. In this section, we first introduce the 
Hamiltonian and the modified second quantization, then derive the heat current density 
operator which includes both the normal velocity and the anomalous velocity from the 
Berry-phase effect. Using the Green-Kubo formula, the general formula of the phonon 
Hall conductivity is obtained. Then we give an application example on the kagome 
lattice in Sec. El In this section the computation details about the dynamic matrix, 
the Chern numbers and the phonon Hall conductivity are given, and the behaviors and 
relations between the phonon Hall conductivity, Chern numbers, and the band structures 
are discussed. In the end a short conclusion is presented in Sec. HI 

2. The PHE theory 

In this section, we give the detailed derivation for the theory of the PHE. We use the 
Hamiltonian in Refs. [26] and [31], which is a positive definite Hamiltonian to describe 
the ionic crystal lattice with in an applied magnetic field. 

2.1. The Hamiltonian and the second quantization 

The Hamiltonian for an ionic crystal lattice in a uniform external magnetic field 
|261 l28l [31] can be written in a compact form as 

I ~ ~ 1 

H = -{p- Auf{p - Au) + -u^Ku 

II _ _ 

= -p^p+ -u^{K - A^)u + u^Ap. (1) 
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Here, m is a column vector of displacements from lattice equilibrium positions for all the 
degrees of freedom, multiplied by the square root of mass; p is the conjugate momentum 
vector, and K is the force constant matrix. The superscript T stands for the matrix 
transpose. A is an antisymmetric real matrix, which is block diagonal with elements 

A = ^ II 1^ ^ (in two dimensions), where h is proportional to the magnetic field, 

and has the dimension of frequency. For simplicity, we will call h magnetic field later. 
According to [9], h is estimated to be 0.1 cm~^ 3 x lO^radHz at a magnetic field 
B = n and a temperature T = 5.45 K, which is within the possible range of the 
coupling strength in ionic insulators [SI |35] . The on-site term, u^Ap, can be interpreted 
as the Raman (or spin-phonon) interaction. Based on quantum theory and symmetry 
consideration, the phenomenological description of the spin-phonon interaction was 
proposed many years ago [3ll [351 [321 [Ml [361 [33 1311 [31]. From the first row of Eq. ([T|), 
we find both of the two terms are positive definite, thus the Hamiltonian ([T]) is positive 
definite. The origin of the Hamiltonian for the PHE is discussed in detail in the 
supplementary information of Ref. [26] . 

The Hamiltonian Eq. ([T]) is quadratic in u and p. We can write the linear equation 
of motion as 

p= -{K~ A'^)u - Ap, (2) 

u = p — An. (3) 

The equation of motion for the coordinate is, 

u + 2Au + Ku = Q. (4) 

Since the lattice is periodic, we can apply the Bloch's theorem ui = ^e'^i'^r^-^'t) _ rjj^^ 
polarization vector e satisfies 

[{-iio + Af + D]e = 0, (5) 

where D(k) = —A? + Knie^'^^'"'^^^'^ denotes the dynamic matrix and A is block 
diagonal with elements A. D^Ki^y^ and A are all nd x nd matrices, where n is the 
number of particles in one unit cell and d is the dimension of the vibration. 
From Eq. ([5|), we can require the following relations: 

elfc = efc; uj-k = -Wfc- (6) 

Here, we use the short-hand notation k = (k, cr) to specify both the wavevector and 
the phonon branch, and —k means (— k, — a). In normal lattice dynamic treatment, 
we usually take a, w > as a convention, and require cUa^k = oOa-k- For the current 
problem, this is not true [261 [31] ■ It is more convenient to have the frequency taking 
both positive and negative values and require the above equation ([6|). And from Eq. ([3]), 
the momentum and displacement polarization vectors are related through 

Hk = -iuJkf^k + Aek- (7) 

Equation ([5|) is not a standard eigenvalue problem. However, we can describe the 
system by the polarization vector x = (/i, e)"^, where fi and e are associated with the 
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momenta and coordinates, respectively. Using Bloch's theorem, Eqs. ([2]) and (E]) can be 
recasted as: 

Here the Ind is the nd x nd identity matrix. Therefore, the eigenvalue problem of the 
equation of motion ([8]) reads: 

■^k '^k •) '^k HeS = Wfc Xfc . (9) 

where the right eigenvector Xk = (/ifc, Cfc)^, the left eigenvector xj = (e|., —fJ^],)/ {—2iujk)-, 
in such choice the second quantization of the Hamiltonian Eq. ([1]) holds, which will be 
proved later. Because the effective Hamiltonian H^s is not hermitian, the orthonormal 
condition then holds between the left and right eigenvectors, as 

5a,k ^^T'.k = 5a<7'- (10) 

We also have the completeness relation as 

X! ^<T,k ® 5 J,k = ^2nd- (11) 

cr 

The normalization of the eigenmodes is equivalent to [11] 

eiek + —e\Aeu = l. (12) 

From the eigenvalue problem Eq. ([9]), we know that the completed set contains the 
branch of the negative frequency. And from the topological nature of the PHE [26], 
the formula of the phonon Hall conductivity can be written in the form comprises the 
contribution of all the branches including both positive and negative frequency branches. 
In order to simplify the notation, for all the branches, we define 

a-k = 4- (13) 
The time dependence of the operators is given by: 

ak{t) = aue-'^^\ (14) 

4(t) = 4e^"'=*. (15) 
The commutation relation is 

[afc,4'] = 4,fc'sign(o-). (16) 
And we can get 

(4«fc) = /(t^fc)sign(cT); (17) 
(afc4) = [l + /(u;fc)]sign(a). (18) 
Here f{ujk) = (e^'^'^/'^'^-sT) _ -[^-i jg ^]-^g Bose distribution function. 
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The displacement and momentum operators can be written in the following second 
quantization forms 



(19) 



Here, \u:k\ = ujksign{a). We can verify that the canonical commutation relations are 
satisfied: [ui,p]/] = ih6u'Ind by using the completeness Eq. (fTT!) and the commutation 
relation Eq. ( fT6l) . The Hamiltonian Eq. ([T]) then can be written as [H] 

1 / ASi,u Ki^v - A^di^v \ 

H=7^2^Xl \ 'jr 'jr ]Xl' (21) 



where 



^'=\Z,) =\\^.Y.^^^""'-''^^^^-^ (22) 
X.= (_^;) =^^E5^^e--'-5^,4. (23) 



Here = y^j^ and = {-2iujk)^^^y It is easy to verify that \xu Xr] = -i^^whnd- 
Because of e*^^'' '''"^' '^) = and the definition of the dynamic 

matrix D, then the Hamiltonian can be written as 

k,k',i 

= ^II^I^^fcl4«fc> (24) 

k 

which contains both the positive and negative branches. Here we use the identity 
g«R-r(k'-k) _ and the eigenvalue problem Eq. ([9]). Using the relations 

Eqs. (IT51) and (ITBi) . it is easy to prove that Eq. is equivalent to the form 

H = J2a>o,k ^^kicik^k + 1/2) which only includes the nonnegative branches. 

2.2. The heat current operator 

The heat current density can be computed as [IQ] : 

J = :^T.i'^i-^i'WKii'Ui', (25) 
'^^ 1,1' 

where V is the total volume of unit cells. Because of the equation of motion Eq. ([3]), 
we can rewrite the heat current as 

3 = ^T.xlMn'Xi', (26) 

1,1' 
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Inserting the Eqs. fl22|23l) . we obtain 



Because of 



^ = Wn ^ CkCk'e'^^''-'''-^^-''^xlMii,Xk'aiak' 

k,k'Ll' 



^ ^.R,(ic'-ic) j2 e^(^.'-i^^)-'^'(R, - R,)Ka> = ^N6^>^^, (29) 

I I' 



the heat current can be written as 

if}/ — ^ ^ ^rp dH -j- { QC\\ 

^ ^ Cf7,kCCT',ka:;<T,k a;<j',k«'^ka<^',k, (30) 



(T,(T',k 

here we use 



I f 



(31) 



by making the first derivative of the square of the effective Hamiltonian Eq. ([H]) with 
respect to the wave vector k. From the eigenvalue problem Eq. ([9]), we have 

H^^X = Xn; X^H^s = ^X^- (32) 

Where the 2nd x 2nd matrices X = {xi,X2, ...,X2nd) = {^a} (the system has 2nd 
phonon branches), X = {x^}, and Q = dia.g{uji,uj2, ...,uj2nd) = Wa}- Because of the 
completeness relation Eq. ( ITTl) . XX'^ = l2nd, we get 

H^s = Xn^X^. (33) 

By calculating the derivative of the above equation, and using the definition of Berry 
connection, 

.4^A-f. (34) 
Taking the first derivative of Eq. (133|) with respect to k, we obtain 

S^.x{^HA,^^i)x^. (35) 

Because of the orthogonality relation between left and right eigenvector Eq. ( ITOl) . at last 
we obtain the heat current as 

ih ( dVl? \ 

(T,(T',k \ / (T,(j' 

The first term in the bracket is a diagonal one corresponding to Wo-^f relating the 
group velocity. The second term in the bracket [A, Vt^] gives the off-diagonal elements 
of the heat current density, which can be regarded as the contribution from anomalous 
velocities similar to the one in the intrinsic anomalous Hall effect. The Berry connection 
A, or we can call it Berry vector potential matrix (the Berry vector potential defined 
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in Ref. A'^(k), is equal to iA'^'^ = ix^^-), induces the anomalous velocities to the 
heat current, which will take the responsibility of the PHE. Therefore, the Berry vector 
potential comes naturally into the heat current and the PHE. Such a picture is clearer 
than that in Ref. I26l. 



2.3. The phonon Hall conductivity 

Inserting the coefficients c and c to Eq. ( |36ll . we get 



This expression is equivalent to that given in Refs. |26] and |3T]. Based on such 
expression of heat current, the phonon Hall conductivity can be obtained through the 
Green-Kubo formula 1^: 



y rn/{kBT) 

K^y = —j^ dX dt {r{-^X)Jy{t))^^, (38) 

where the average is taken over the equilibrium ensemble with Hamiltonian H. The 
time dependence of the creation and annihilation operators are given as Eqs. ( ITil) and 
(fTSjl . which are also true if t is imaginary. From the Wick theorem, we have 

(«I,k«<T',ka^a^',k) = («I,k'^<T',k)(a^a^/ k) 

+ (aU4,k)(^'^'.k«a',k) 

+ (4,ka<f',k)(aa',k4,k)- (39) 

Using the properties of the operators and a as Eq. (ITTl) . we have 
(4,k«^',k)(4,k^^"',k) 

= f{uja,k)fi^a,^)Saa'S^-,sign{a)sign{a), 
(4,kat k)(aa',ka^-',k) 

= /(wa,k)(/(wa',k) + I)5k,_k5<7,-a5a' -a'Sign(a)sign(a') 

(4,k«<T"',k)(aa',k4,k> 

= /(wa,k)(/(c^a',k) + l)(5k,k'^^v^<T',<7sign(a)sign(a'). 

Similar as that in Ref. [2B], the diagonal term in the bracket corresponding to 
Wo-^^ has no contribution to the phonon Hall conductivity because which is an odd 
function of k. Because of the off-diagonal term 

[A,^,nX,^, = iul,-ul)AZ' (41) 
and A^'^' = x'^^r- from the definition, the phonon Hall conductivity can be written as 

h 

k,(J,(T'^(T 

. pt t 



(40) 
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Here we simplify the notation of the subscripts of e which have the same wave vector 
k. We can prove K^y = —Ky^, such that 

h 



«.„ = r^ E [/K)-/K')](^<x + c^.')Xt' (43) 



here 



l- 



= - t {Ai:'Al'; - {k, ^ ky)) , (44) 

in the last step we use the relation x^^^Xo-' = (wo-' — ^a-)x^^Xa' and the definition 
of A in Eq. ( I34p . And the Berry curvature is 

Rcr \ ^ RCCr' 

-^h h / -^h h 



'^Y.{^lxK-ikx^ky)) 

cr' 



(45) 



The definition of Berry curvature here is the same as that of Ref. 126], that is, 
BJ,,, = J_AJ. - J-aI Rom the above der.vatio,, we flud that a BerT^ curvature 

can be defined uniquely for each band by looking at the phases of the polarized vectors 
of both the displacements and conjugate momenta as functions of the wave vector. If 
we only look at the polarized vector e of the displacement, a Berry curvature cannot 
properly be defined. We need both e and /i. The nontrivial Berry vector potential takes 
the responsibility of the PHE. The associated topological Chern number is obtained 
through integrating the Berry curvature over the first Brillouin zone as 

where, L is the length of the sample. 



3. Application on the kagome lattice 

In Ref. [23, we provide a topological understanding of the PHE in dielectrics with 
Raman spin-phonon coupling for the honeycomb lattice structure. Because of the nature 
of phonons, the phonon Hall conductivity, which is not directly proportional to the Chern 
number, is not quantized. We observed a phase transition in the PHE, which corresponds 
to the sudden change of band topology, characterized by the altering of integer Chern 
numbers. Such PHE can be explained by touching and splitting of phonon bands. To 
check whether the mechanism of the PHE is universal, in the following we apply the 
theory to the kagome lattice, which has been used to model many real materials |30] . 
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Figure 1. (color online) The schematic picture of kagome lattice. Each unit cell has 
three atoms such as the number shown 1,2,3. The coupling between the atoms are 
Km, Ko2, Ko3. Each unit cell has six nearest neighbors; the coupling between the unit 
cell and the neighbors are Ki, K2, Kq. 



3.1. Calculation of the dynamic matrix D 

In order to calculate the phonon Hall conductivity, we first need to calculate the dynamic 
matrix -D(k), for the two-dimensional kagome lattice. As shown in Fig. [H each unit cell 
has three atoms, thus n = 3. We only consider the nearest neighbor interaction. The 
spring constant matrix along x direction is assumed as 

Kl = 0.144 eV/(uA^) is the longitudinal spring constant and the transverse one Kt is 
4 times smaller. The unit cell lattice vectors are (a, 0) and (a/2, 0^3/2) with a = lA. 

To obtain the explicit formula for the dynamic matrix, we first define a rotation 
operator in two dimensions as: 



Ui9) 



cos 6 — sin 6 
sin 6 cos 6 



The three kinds of spring-constant matrices between two atoms are Kqi = (between 
atoms 1 and 2 in Fig. [1]), Kq2 = t/(7r/3)i^^f/(— 7r/3) (between atoms 2 and 3), 
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Figure 2. (color online) The contour map of dispersion relations for the positive 
frequency bands. For all the insets, the horizontal and vertical axes correspond to 
wave vector and ky, respectively. The upper six insets are the dispersion relations 
for bands 1 to 6 (from left to right) at ft = 0, respectively. And h = 10 rad/ps for the 
lower ones. 



= U{—TT/3)KrcU{7i/3) (between atoms 3 and 1), which are 2x2 matrices. Then 
we can obtain the on-site spring- const ant matrix and the six spring-constant matrices 
between the unit cell and its nearest neighbors as: 



/ 



V 



2(i^^oi + ^02) 
-K02 



2iKoi + K03] 

-Ko3 



2{Ko2 + Ko3) J 





-^01 




\ 

^ 

/ \ 


-Ko3 

/ -7^02 \ 








( -J^oi 



\ 

^ 




^ 

/ \ 

-7^03 





which are 6x6 matrices. Finally we can obtain the 6x6 dynamic matrix -D(k) as 

D(k) = -A^ + Kq + Kie'^^ + K2e'^-^ ' 
+ K3e'^-^+^^ + K^e-'^^ 



(48) 



where, = —h'^ ■ Iq, here Iq is the 6x6 identity matrix. 
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3.2. The PHE and the associated phase transition 

After we get the expression for the dynamic matrix, we can calculate the eigenvalues and 
eigenvectors of the effective Hamiltonian. Inserting the eigenvalues, eigenvectors and the 
D matrix to the formula Eq. (H3|) , we are able to compute the phonon Hall conductivity. 
As is well known, in quantum Hall effect for electrons, the Hall conductivity is just 
the Chern number in units of e^/h {h is the Planck constant); thus with the varying 
of magnetic field, the abrupt change of Chern numbers directly induces the obvious 
discontinuity of the Hall conductivity. However, for the PHE, there is an extra weight 
of {uj(j + Wo-')^ in Eq. ( 143|) . which can not be moved out from the summation. As 
a consequence, the change of phonon Hall conductivity is smoothened at the critical 
magnetic field. However, in the study on the PHE in the honeycomb lattice system 
[26], from the first derivative of phonon Hall conductivity with respect to the magnetic 
field h, at the critical point he, we still can observe the divergence (singularity) of 
dKxy/dh, where the phase transition occurs corresponding to the sudden change of the 
Chern numbers. Can such mechanism be applied for the kagome lattice system? In the 
following, we give a detailed discussion on it. 

Inserting the dynamic matrix Eq. ( I48l) to the effective Hamiltonian Eq. (jO]), we 
calculate eigenvalues and eigenvectors of the system, and also get the dispersion relation 
of the system. Because each unit cell has three atoms, and we only consider the two- 
dimensional motion, we get six phonon branches with positive frequencies. The branches 
with negative frequencies have similar behavior because of U-k = —^k- We show the 
contour map of the dispersion relation in Fig. O We can see that the dispersion relations 
have a 6-fold symmetry. For different bands, they are different. With a changing 
magnetic field, the dispersion relations vary. The point F (k = (0,0)) is the 6-fold 
symmetric center; the point K (k = (^,0)) is 3-fold symmetric center; and the middle 
point of the line between two 6- fold symmetric centers, X (k = (vr, ^^^)) is a 2- fold 
symmetric center. In the following discussion, we will see the possible bands touching 
at these symmetric centers. 

Using the formula Eq. ( H3|) . we calculate the phonon Hall conductivity of the kagome 
lattice systems, the results are shown in Fig. [31 Similar as shown in Ref . [26] , we find the 
nontrivial behavior of the phonon Hall conductivity as a function of the magnetic field. 
When h is small, Hi^y is proportional to h, which is shown in the inset of FigjSl while the 
dependence becomes nonlinear when h is large. As h is further increased, the magnitude 
of Kxy increases before it reaches a maximum magnitude at certain value of h. Then the 
magnitude of K^y decreases and goes to zero at very large h. The on-site term in the 
Hamiltonian ([T]) increases with h quadratically so as to blockade the phonon transport, 
which competes with the spin-phonon interaction. Because of the coefficient of f{uJa) 
in the summation of the formula Eq. (H3|) . the sign of the Hall conductivity will change 
with temperatures, which is clearly shown in the inset of Figl3l While the phonon hall 
conductivity at weak magnetic field is always positive for the honeycomb lattice, the 
sign reverse of the phonon Hall conductivity with temperature for the kagome lattices 
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Figure 3. (color online) The phonon Hall conductivity vs magnetic field at different 
temperatures. The inset is the zoom-in curve of the phonon Hall conductivity at weak 
magnetic field. Here the sample size iVi=400. 



is novel and interesting, which could be verified by future experimental measurements. 

We plot the curves of the Chern numbers of bands 2 and 3 as a function of the 
magnetic field in Fig. IH The phonon Hall conductivity at T = 50K is also shown 
for comparison. To calculate the integer Chern numbers, large number of k-sampling 
points is needed. However there is always a zero eigenvalue at the F point of the 
dispersion relation, which corresponds to a singularity of the Berry curvature. Therefore, 
we cannot sum up the Berry curvature very near this point to obtain Chern number 
of this band, unless we add a negligible on-site potential |M^V"onsiteW to the original 
Hamiltonian |26], which will not change the topology of the space of the eigenvectors. 
In Fig. m we set V^nsite = 10~^Kl. The Chern numbers of bands 2 and 3 have three 
jumps with the increasing of the magnetic field, although the phonon Hall conductivity 
is continuous. For other bands, the Chern numbers keep constant: = = — 1, 
= 0, and = 1. For the electronic Hall effect, we know it is quantized because 
the Hall conductivity is directly proportional to the quantized Chern numbers. Here 
we also find the quantized effect of the Chern numbers from Fig. [3l while there is no 
quantized effect for the phonon Hall conductivity. Such difference of the PHE from the 
electronic Hall effect comes from the different nature of the phonons respective to the 
electrons. In Eq. (H3l) . in the summation, an extra term [ua + cUo-')^ relating to the 
phonon energy which is an analog of the electrical charge term in the electron Hall 
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Figure 4. (color online) The Chern numbers and the phonon Hall conductivity vs 
magnetic field. The dashed line and the dotted line correspond to the Chern numbers 
of phonon bands 2 and 3 (left scale). The solid line correspond to the phonon Hall 
conductivity (right scale) at T = 50 K. 

effect, can not be moved out from tlie summation. Combining tlie Bose distribution, 
tlie term f{ujcr){u)a + Wo-/)^ make tlie plionon Hall conductivity smooth, no discontinuity 
comes out although the Chern numbers have some sudden jumps. From the discussion 
in Ref . [26] , the discontinuity of the Chern numbers corresponds to the phase transitions 
and would relate to the divergency of derivative of the phonon Hall conductivity. 

Figure [5] shows the curves of the derivative of the phonon Hall conductivity and 
the Chern numbers at the critical magnetic fields. The first derivative of phonon Hall 
conductivity has a minimum or maximum at the magnetic fields hd = 5.07, hc2 = 
6.75,and/ic3 = 20.39 rad/ps for the finite-size sample (the sample has = Nl unit 
cells). The first derivative dn^y/dh at the points h^, hc2, hcs diverges when the system 
size increases to infinity [26] • At the three critical points the second derivative d'^n^y/dh'^ 
is discontinuous, which is shown in the inset of Fig. |5l across which phase transitions 
occur. For different temperatures, the phase transitions occur at exactly the same 
critical values. Thus the temperature-independent phase transition does not come from 
the thermodynamic effect, but is induced by the topology of the phonon band structure, 
which corresponds to the sudden change of the Chern numbers. While there is one 
discontinuity of the Chern numbers for the honeycomb lattice system, for the kagome 
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/7(rad/ps) 

Figure 5. (color online) The first derivative of the phonon Hall conductivity dkxy/dh 
at T = 50K and the Chern numbers of bands 2 and 3 in the vicinity of the magnetic 
fields. The solid line correspond to the dkxy/dh atT ~ 50K (left scale); the dashed and 
dotted lines correspond to the Chern numbers of bands 2 and 3, respectively (right 
scale). The inset shows the second derivative with respective to the magnetic field 
dk^y/dh^ (vertical axis) vs magnetic field h (horizontal axis) at T = 50 K. 

lattice system, there are three ones corresponding to the divergency of the derivative of 
the phonon conductivity, which can be seen in Fig. [51 

The touching and sphtting of the phonon bands near the critical magnetic field 
induces the abrupt change of Chern numbers of the phonon band p6]. In Ref. |26j . 
for the PHE in the honeycomb lattices, we know that band 2 and 3 are going to touch 
with each other at the F point if the magnetic field increases to h^, at the critical 
magnetic field, the degeneracy occurs and the two bands possess the cone shape; above 
the critical point h^ the two bands split up. Therefore, the difference between the two 
bands decreases below and increases above the critical magnetic field, and is zero at the 
critical point. The eigenfrequecy difference is in the denominator of the Berry curvature, 
thus the variation of the difference around the critical magnetic field, dramatically affects 
the Berry curvature of the corresponding bands. In the kagome lattice systems, we find 
that the touching and splitting of the phonon bands not only occurs at the T point, but 
also occurs at other points, which is shown in Fig. [6l At the first critical points hd, the 




Figure 6. (color online) The dispersion relations around the critical magnetic fields, 
(a), (b), and (c) show the dispersion relations along the direction from F (k=(0,0)) to K 
(k ~ (^,0)) and to X (k = (tt, -^^)) at the critical magnetic fields ft-ci = 5.07rad/ps, 
hc2 = 6.75rad/ps, and /ics = 20.39rad/ps, respectively, (d)-(f) show the contour maps 
of the dispersion relation of band 2 at the three critical magnetic fields, (g)-(i) show 
the contour maps of the dispersion relation of band 3 at the three critical magnetic 
fields. The squares with number 1, 2, and 3 are marked for the touching points. In (d), 
(g) and (e), (h), we only mark one of the six symmetric points by squares of number 
1 and 2 for simplicity. 
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bands 2 and 3 touch at the point K (marked by a square with number 1); at hc2 the two 
bands touch at X ( marked by a square with number 2); while only for the third critical 
one hcZi band 2 and 3 degenerate at the point T (marked by a square with number 3). 
From the contour maps of bands 2 and 3, we clearly see that the critical magnetic fields 
hcii hc2, and hcs, there are local maximum for the band 2 and the local minimum for the 
band 3. Thus for all the critical magnetic fields where the Chern numbers have abrupt 
changes, in the wave-vector space we can always find the phonon bands touching and 
splitting at some symmetric center points. 

Therefore, through the study of the PHE in both honeycomb lattices [26] and 
kagome lattices, we find discontinuous jumps in Chern numbers, which manifest 
themselves as singularities of the first derivative of the phonon Hall conductivity with 
respect to the magnetic field. Such associated phase transition is connected with the 
crossing of band 2 and band 3, which corresponds to the touching between a acoustic 
band and a optical band. However, we can not observe the similar associated phase 
transition in triangular lattices because of no optical bands, where the Chern number 
of each band keeps zero while the phonon Hall conductivity is nonzero because of the 
nonzero Berry curvatures. 

4. Conclusion 

We present a new systematic theory of the PHE in the ballistic crystal lattice system, 
and give an example application of the PHE in the kagome lattice which is a model 
structure of the many real magnetic materials. By the proper second quantization for the 
Hamiltonian, we obtain the formula for the heat current density, which considers all the 
phonon bands including both positive and negative frequencies. The heat current density 
can be divided into two parts, one is the diagonal, another is off-diagonal. The diagonal 
part corresponds to the normal velocity; and the off-diagonal part corresponds to the 
anomalous velocity which is induced by the Berry vector potential. Such anomalous 
velocity induces the PHE in the crystal lattice. Based on such heat current density 
we derive the formula of the phonon Hall conductivity which is in terms of the Berry 
curvatures. From the application on the kagome lattices, we find that at weak magnetic 
field, the phonon Hall conductivity changes sign with varying temperatures. It is also 
found that the mechanism on the PHE about the relation between the phonon Hall 
conductivity, Chern numbers and the phonon band structure can be generally appUied 
for kagome lattices. While there is only one discontinuity in PHE of the honeycomb 
lattices, in the kagome lattices there are three singularities induced by the abrupt change 
of the phonon band topology, which correspond to the touching and splitting at three 
different symmetric center points in the wave-vector space. 
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